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ABSTRACT 

We present multidimensional non-LTE radiative transfer models of hydrogen and helium 
line profiles formed in the accretion flows and the outflows near the star-disk interaction re- 
gions of classical T Tauri stars (CTTSs). The statistical equilibrium calculations, performed 
under the assumption of the Sobolev approximation using the radiative transfer code TORUS, 
has been improved to include He I and He II energy levels. This allows us to probe the physical 
conditions of the inner wind of CTTSs by simultaneously modelling the robust wind diagnos- 
tic line He I A10830 and the accretion diagnostic lines such as Pa/3, Br7 and He I A5876. The 
code has been tested in one and two dimensional problems, and we have shown that the results 
are in agreement with established codes. We apply the model to the complex flow geometries 
of CTTSs. Example model profiles are computed using the combinations of (1) magneto- 
spheric accretion and disc wind, and (2) magnetospheric accretion and the stellar wind. In 
both cases, the model produces line profiles which are qualitatively similar to those found in 
observations. Our models are consistent with the scenario in which the narrow blueshifted 
absorption component of He I A10830 seen in observations is caused by a disc wind, and the 
wider blueshifted absorption component (the P-Cygni profile) is caused by a bipolar stellar 
wind. However, we do not have a strong constraint on the relative importance of the wind and 
the magnetosphere for the 'emission' component. Our preliminary calculations suggest that 
the temperature of the disc wind and stellar winds cannot be much higher than ^10,000 K, on 
the basis of the strengths of hydrogen lines. Similarly the temperature of the magnetospheric 
accretion cannot be much higher than 10,000 K. With these low temperatures, we find that 
the photoionzation by high energy photons (e.g. X-rays) is necessary to produce He I A10830 
in emission and to produce the blueshifted absorption components. 
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1 INTRODUCTION 

Mass-loss processes in young stellar objects (YSOs), such as the 
classical T Tauri stars (CTTSs), are a fundamental problem in the 
star formation process. The origin and the physical conditions of 
the wind from YSOs are crucial to our understanding of the proto- 
stellar phase and of the evolution of protostellar discs. The mass- 
loss process is closely related to the evolution of the stellar rotation 
since it is one of the possible mechanisms through whi ch YSOs 
may lose their angular momentum (e.g . iHartmann & Stauffed 
1 19891 : iMatt & Pudiit3l2005l : l2007j : l2008iJ) . The mass-loss due to 
winds may solve the so-called 'angular-momentum problem', in 
whic h a large fraction o f accreting CTTS have rather slow rota- 
tions jHerbst et al]|2007h at ~ 10 per cent of break-up speed (e.g., 
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iMatt & Pudrit j|2007lj) despite of a rather large amount of the an- 
gular momentum added to stars by accretion process. 

Most likely, the outflows are formed in magnetohydro- 
dynamics processes in which a large-scale open magnetic 
field is anchored to a rotating object; however, whether 
the o bject is the disk, central star or both is still un- 



clear Edwards_et^l. 2006; JFeroira^ I2OO6I ; 
iKwan. Edwards. & Fischeil l2007h . There are at least three pos- 
sible configurations of the wind formation: (1) a disc wind in 
which one assumes a sufficient magnetic field and ionization 
fraction to launch magneto- centrifugal winds over a relatively 
wide range of disk radii, f rom an inner truncation radius out to 
a several au (e. g. Ustvug ova et al.l 19951 ; iRomanova et al. 199^; 



Ouved & Pudrit z.il997i; lUs tvugova e t al.Ml999i; ,Konigl & Pudritzl 
2000; Krasnopol skv et alj2 003; Pudri tz et alj2007h . (2) an X-wind 
in which the wind launching region is restricted to a narrow re- 
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gion near the inner edge of the accretion disk around the coro- 
tation ra dius where the d isk interfaces with the stellar magne- 
tosphere ( IShuetal.lll994h . and (3) a stellar wind in which out- 



face fe.s., Hartmann, Avrett, & Edwards 1982: Kwan & Tademaru 


Il988l; iHirose et alj|l997,; .Romanova et aL 


,20051; iCranmer 


i2009D. 
J2009b 


Interestingly, recent MHD simulations by 


Romanova et al 



found yet another type outflow called 'conical winds' which are 
produced when the stellar dipole magnetic field is squeezed by 
the disc into the X-wind type configuration. The resulting out- 
flow occurs in rather narrow conical shapes with their half open- 
ing angles between 30° to 40° . This conical wind is similar to the 
X-wind mentioned above, but the former is driven by the mag- 
netic force and it does not required to have the magnetospheric 
radius matched with the corotation radius of the system. More re- 
cent simulations performed with larger radial domains have shown 
that the conical winds becom e strongly coUimated by the mag- 
netic force at larger dista nces (Lii. Romanova. & Lovelacg|201ll : 
iKonigl. Romanova. & Lovelace, 20 1i[ k 

In order to distinguish between the different theoretical 
models, and to understand the origin of the angular momen- 
tum transfer, constraints on the wind launching location are 
very important. Spectroscopic studies of strong emission lines 
in CTTSs may help to so l ve this problem. Th e observations 
by e.g . [jakami et ai l j2002h. lEdwards et all ( |2003|) , lOupree et all 
(/005) and Edwards et al.' ('2006') demonstrated a robustness of op- 
tically thick He I A10830 as a diagnostic of the inner wind from 
the accreting stars. Unlike other strong emission lines (e.g. Ha and 
H/?) which also often show a sign of the inner winds. He I A10830 
from CTTS shows P-Cygni type profile (with deep blueshifted 
absorption below continuum level) which resembles that of the 
hot s te llar winds expanding in radial direction (e.g., iKwan et al] 
l2007b . lEdwards et"S] ( l2006h showed that about 70 per cent of 
CTTS exhibit a blueshifted absorption (below continuum) in con- 
trast to Hof which shows about only 10 per cent of stars show simi- 
lar ty p e of absor ption component (e.g..'Reipurt h. Pedrosa. & Lagol 
Il996h . lEdwardsl ( 12007.) and .Kwan et al. ^2007^ suggested that the 
blueshifted absorption component in the He I A10830 profiles is 
cased by a stellar wind in about 40 per cent, and by a disc wind in 
about 30 per cent of the sample s in Edwards et al. (2006). Recent 
local excitation calculations of lKwan&Fischej ( l201ll) in combi- 
nation with the spectroscopic observations demonstrated the effec- 
tiveness of He I A10830 for probing the density and temperature of 
the inner wind. 

In our previous radi at ive transfer models, dHarries 



2005 



20001 : iKurosawa et"ai] |2004 [Svm in gton. Harries. & Kurosawa 



Kurosawa, Harries, & Sv mingtonl 20061 : 
Kurosawa. Rom anova, & Harries 2008), we have considered 



the statistical equilibrium of hydrogen atoms only. In light of the 
recent re cognition of He I A10830 line as a useful wind diagnostic 
tool ( e.g. lEdwards et alll2006l : iKwan et alJboOTi : iKwan & FischeJ 
I20T l|), we improve our model to include helium ions (He I, 
Hell and He III). This will allow as to investigate the physical 
conditions (e.g. geometry, temperature, density and kinematics) 
of the innermost part of the wind more effectiv ely, and will be 
a sign ificant improvement to the wind study of iKurosawa et al.l 
ilQQ& i who used only Ha. If we are to reliably determine mass loss 
rates and mass accretion rates of CTTSs, a more self-consistent 
modelling approach is required. This includes a multiple line fitting 
of observations including both hydrogen and helium to model the 
wind and magnetospheric accretion flows simultaneously. 

Our main goals in this paper are to present the methods and 



tests for our improved radiative transfer code, and to demonstrate 
the capability of the code in the applications to the wind and mag- 
netospheric accretion problems of CTTSs. Detailed investigations 
of the wind and accretions flows, e.g. constraining flow geometry, 
density, temperature and relative strength of the mass-loss to mass- 
accretion rates, will be presented in a future paper as these require 
a large parameter space survey. 

In Section |2l the model assumption, the method of statisti- 
cal equilibrium and the observed profile calculations are presented. 
The tests of the code in one and two dimensional (ID and 2D re- 
spectively) geometries will be given in Section|3] Simple kinematic 
models of a disc wind and a stellar wind in combinations with a 
magnetospheric accretion model will be presented in Section|4| Ex- 
ample profiles of hydrogen and helium are also given in the same 
section. Finally main findings and conclusions are given in Sec- 
tion|6| 



2 RADIATIVE TRANSFER CODE 



The radiative transfer code TORUS ( lHarries |l2000l: 'Kurosawa et al.l 

2004 \ [Svmington et all I2OO5I : [Ku rosawa et al. 200a 

Kurosawa et al. 2008, but see also Acreman, Harries, & Rundlg 
,201Qi : ,Rundle et al..,201Q) is extended to include helium ions (He I, 
Hell and Hem). Our model uses the three-dimensional (3D) 
adaptive mesh refinement (AMR) grid, and it allows us an accurate 
mapping of an MHD simulation data on to the radiative transfer 
grid (without a loss of resolution). The code also works in 2D 
(axisymmetric) and ID (spherically symmetric) cases. The basic 
steps for computing the line profiles are as follows: (1) mapping 
of the density, velocity and temperature, from either an MHD 
simulation output or an analytical model, to the radiative transfer 
grid, (2) the source function (Si,) calculation and (3) the observed 
flux/profile c alculation. In the seco nd step, we adopt the method 
descr i bed bvlKlein & Casto3 (Il978 ) (see also iRvbicki & Hummed 
1 19781 : iHartmann. Hewett. & CalveJ Il994h in which the Sobolev 
approximation method is used. The Sobolev approximation, is 
valid when (1) a large velocity gradient is present in the gas flow, 
and (2) the intrinsic line width is negligible compared to the 
Doppler broadening of the line. In the following, we describe each 
step of the line profile computation in more detail. 



2.1 Grid construction 

Our radiative transfer code has been developed to handle length 
scales in many order of magnitudes and mul ti-dimensional prob- 



lems with out a particular symmetry (e.g., [K urosawa & Hillier 
2001: Svm ington et al.ll2005l : iHarries et alTl200 4: Kurosawa et aL 
20041) . When the gradient or dynamical range of physical quanti- 
ties such as density is very large, a logarithmic scale in the radial 
direction could be used to increase an accuracy of optical depth cal- 
culations if the system is spherically or axially symmetric. When a 
problem requires an arbitrary geometry, there is no simple way to 
construct a logarithmic scale. However, an efficient cubic or square 
grid (depending on the problem is 3D or 2D) can be constructed 
by subdividing a cube/square into smaller cubes/squares where the 
value of density or opacity is larger than a given threshold value. 
We use the 8-way/4-way tree (octal/quad tree) data structure to con- 
struct the grids and to store the physical quantities required in the 
radiative transfer calculations. The data structure used here is sim- 
ilar to a binary tree in which one node splits into two child nodes, 
but ours splits into eight (for 3D) or four (for 2D). 
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for hydrogen are same as in iHarriesI fcOOOl) and ISvmington et al 



Figure 1. An example of the 2D AMR grid constracted for the density from 
the magnetospheric accretion model described by Hartma nn et alj jl994) . 
The grid structure is overlaid on the colour map of the input density. The 
density map are logarithmically scaled and in cgs units. The axes are in 
10^" cm. Although not shown here, the accreting star is located at {x, y) = 
(0, 0), and its radius is 1.4 X lO'^^ cm. 



The algorithm us ed to construct the grid in this paper is very 
similar to that used in lKurosawa &Hillieil(l200 l'). Starting from a 
large cubic/square cell (with size D) which contains a whole com- 
putational domain, we first compute the average density of the cell 
by randomly choosing 100-1000 points in the cell and evaluating 
the density at each point. The average density is then multiplied by 
the volume of the cell (D^) to find the total mass (Afcdi) in the 
cell. If the mass is larger than a threshold mass (Mth), which is a 
user defined parameter, this cell is split into 8 or 4 subcells with 
size D/2. If the mass of the cell is less than the threshold, it will 
not be subdivided. The same procedure is applied to all the subcells 
recursively, until all the cells contain a mass less than the threshold 
(Mceii < A/th)- Alternatively, if the opacity of the field is known 
prior to the construction of the grid (which is usually not the case), 
one can use the optical depth across a cell (tccii) and use the condi- 
tion, Tceii < Tth where rth ~ 0-1, as a condition for cell splitting. 
More complex algorithms could be used to optimise the compu- 
tation for a specific problem. Once the data structure is built, the 
access/search for a given cell can be done recursively, resulting a 
faster code. Fig.[T]shows an example of the AMR grid constructed 
for t he axisymmetric (dipol ar) magnetospheric accretion columns 
from lHartmann et all ( Il994l) . 



2.2 Statistical equilibrium 

Our statistical equil ibriu m computat ion e ssentially follows the 
met hod described bvlKlein & Cast oj h978l) (hereaf ter KC78; see 
also lRvbicki & Hummeifr978l : iHartmann et aljl994l) who uses the 
Sobolev approximation. A main difference between the model of 
KC78 and ours is in the atomic model. Our atomic model consists 
of 20 bound levels for H I and the continuum level for H II, 19 levels 
for He I (up to the principal quantum number n = A level), 10 levels 
for He II and the continuum level for He III. The atomic data used 



1 2005), and those for helium ions are adopted from lHillier & Miller 



1998i") . Here, we have a rather small set of ions and levels; how- 



ever, the method which are to be described next is not limited to 
these ions. 

Since we intend to apply our model to multi-dimensional 
problems (2D and 3D), the main limiting factor of the size of the 
atomic model is computer memory, and CPU time. For example, 
our typical 2D model with CTTS accretion funnels and wind re- 
quires about 2 X 10* grid points (~ 141^) and about 300 Mb of 
RAM; hence adding more complicated atoms and levels will be 
not feasible without a major modification to the code. Note that our 
code is parallelised under Message Passing Interface (MPI), but the 
domain decomposition of the computational grid has not been im- 
plemented because the code is optimised for speed. Adding more 
ions and modification of the code to handle a grid with domain de- 
composition will be done in the future. 

In the original model of KC78, a simpler atomic model was 
used. They assumed that a majority of the helium atoms were ion- 
ized, and excluded He I levels in their calculation. Their assump- 
tion is acceptable since their objective is to model emission lines 
in the hot wind of Of stars which have TcS >30,000 K. The same 
assumption can not be applied to flows with much lower tempera- 
tures such as in the CTTS accretion flows and wind; hence, we have 
added those extra levels. 

In the following, we summarise our computational steps. The 
following equation of statistical equilibrium is applied to each lev- 
els of H I, He I and He II. The total rates for j-th level can be written 
as: 

Rj — Rj + H,^ — Rj = 0, (1) 

where R^ is the net transition rate between j-th level and the lower 
bound levels, and similarly ii^ is the net transition rate between 
j-th level and the upper bound levels. The last two terms, Rf and 
_Rj, are the net recombination and the ionization rates, respectively. 



The rates in equation ([TJ are: 
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R-i 



R^ 
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R'j = Uj 



i<j 
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X exp \~-^]'^^ ~^ n^Cjk 



4tv 



hv 



(2) 



(3) 



(4) 



(5) 



where Aij and Bij are the Einstein coefficients between levels i 
and j, and dj is the collisional rate coefficient. The first and the 
second subscripts in A, B and C correspond to the initial and the 
final state of a transition, respectively. The subscript k denotes the 
continuum state. The LTE and non-LTE populations of i-th level 
are written as n* and Ui, respectively. While aj (u) is the photoion- 
ization cross section of i-th level, Ji, is the mean (angle averaged) 
intensity of the continuum radiation. The symbol Jij refers to the 
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angle and line profile a veraged intensi ty of the line transition be- 
tween levels i and j (cf . iMihalasIl 1 97 Sl) . The electron number den- 
sity and temperature are specified as rie and T. Other symbols have 
usual meaning. Using the Sobolev approximation (the escape prob- 
ability theo ry), Jij , for t he transitions with i < j, can be simplified 
as (see e.g. lCastoilll970l) : 



2/ii/f. 



5i ^ 

gi rii 



(6) 



where Ic.ij is the intensity of the continuum around the line fre- 
quency Vij . The symbols gi and gj denote the degeneracy (statisti- 
cal weight) of the energy levels. Here, we assume that the 'contin- 
uum emissivity and opacity of the gas are veiy small, and the pho- 
tospheric photons are unaffected by continuum processes. In equa- 
tion (|6), we neglect the contribution from the non-local emi ssion 
which are included in the models of iHartmann et aL 1 ( Il994l) and 
[Muzerolle et al. ( 2001). The effect of this omission is very small 
as we will see later in Section [T2l The escape probabilities Pa 
and P c.i-i in equation Q are (c.f. iCastoJl970l[Rvbicki & Hummen 

[Hzi):' 



and 



-dn, 



, ij — 



/ - 



(7) 



(8) 



respectively. In equation ([§}, the integration is performed over the 
solid angle of stellar photosphere (Sl») as seen from the location 
where Pdj is evaluated. The optical depth Tij is defined as 



din 



(9) 



Here, dvn/dln is the derivative of the velocity field projected in 
the direction n, i.e. Vn = v ■ n, and din is the line element in the 
direction of n. The line opacity Xij is: 



rUeC 



giUj 



(10) 



where fij is the oscillator strength. 

Since we assume that the gas is optically thin to the photo- 
spheric continuum radiation, the following simplifications are used 
for the continuum intensity in equation ^ and the mean intensity in 
equations l|4j and (|5}. We use the p l ane-pa rallel model atmosphere 
(the Eddington flux H,y) of lKuruca ( Il979l) . then they become: 



and 



(11) 



(12) 



where W is the geometrical dilution factor W = [1 ~ 
(1 ~ Rl/r'^Y^^]/2 (cf. lMihalas|[T973) . Alternati vely, we can as- 



sume that the photosphere is a blackbody (as in IHartmann et al.l 
ll994lMuzerolle et al.ll200ll) . then we simply have Ic,ij = (T) 
and Jv ~ WBu (T) with an appropriate photospheric temperature 
T. 

The hydrogen particle number conservation gives 



iVHI 



(HI) + n(HII) = 



pX 



{X + 4Y) 



(13) 



where X and Y are the abundances of hydrogen and helium (by 



number), and nip is the proton mass. The symbol A'^ei denotes the 
number of Hi levels (20). Similarly, the helium particle number 
conservation gives 

n, (He I) + ^ n> (He II) + n (He III) 

i=l i=l 

PY 

{X + AY) mp ' 

where the numbers of energy levels used here are A^hcI = 19 and 
A'^HgII = 10. Unless specified otherwise, we use X = 0.9 and 
y = 0.1 for the abundances. Finally, the charge conservation gives 



(14) 



, (HII) + ^ n, (Hell) + 2n (He HI) = ' 



(15) 



Equations lUll, l |13t , J14t and l |15t form a set of algebraic equations 
(A^Hi + A^Hci + A^Hcii + 3 equations) which contain TVhi + A'^hcI + 
A'neii + 3 unknowns. The unknowns here are the level populations 
of H I, He I, He II respectively, and the last 3 unknowns are the 
continuum levels and the electron density, i.e. n (H II) , n (He III) 
and Tie . The set of non-linear algebraic equations will be solved by 
using the Newton-Raphson iteration method. As a starting point on 
the iteration, we use either (1) a LTE population or (2) a set of non- 
UTE solutions from a nearby point which had converged prior to 
the current computation. The level populations normally converge 
within 30 iterations. 

As mentioned earlier, in our models presented here, the gas 
temperature T is given as an input because the exact physical pro- 
cesses that determine the plasma temperature of the accretion flow 
and outflow around CTTS (within a few stellar radii) are not well 
understood. We therefore hold the temperatures fixed and do not 
solve the equation of radiative equilibrium. 



2.3 Observed profile calculations 

The method for the pr ofile calculation is essentially the same as in 
iKurosawa et al Here, we briefly summarise the steps. Once 

the level populations of H and He are obtained (as in Section l2^ . 
the emissivity and opacities for the continuum and line {rjc, Xc rji 
and y;) can be readily obtained (e.g. equation llOt. By following 
iMihalaj h978l) . the total source function (S^) can be written in 
terms of the continuum and line source function which are defined 
as Sc = r]c/Xc and 5*; =r]i/xi, respectively, i.e. 



(j)i,Si + hSc 



(16) 



where h = Xc/xi and 0,^ is a normalized line profile. When 
a line broadening is negligible (f>v is a norm al Doppler profile - 
Ho wever, as noted a nd de monstrated by Muz erolle et al. I l l200lh 
and iKurosawa et"ai] ( |2006|) . with a moderate mass-accretion rate 
( 10~^ Mq yr~^) of a CTTS, Stark broadening becomes important 
in the optically thick lines, e.g. in Ha and possibly in He I A10830. 
For those lines, we adopt a Voigt profile to include the broad- 
ening effect. The normalized Voigt profile is written as </)iy = 
^-1/2 ^Q,^ y-^ where 



H{a,y) = - 



(y - y') + 



dy' 



(17) 



In the expression above, a — r/47rAi/D, y = {v — i^o) /A/^d, and 
y' = [v' — vq) /Ai^D (c . f . iMihalasI 1 97 st) . vo is the line centre fre- 
quency, and Ai/D is the Doppler line width of an atom (due to its 
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thermal motion) which is given 

where mi is either the mass of a hydrogen atom (ttih) or that 
of a helium atom (mnc). The damping constant F, which de- 
pends on the physical condit ion of the gas, is parametrised by 
iLuttermoser & Johnsonid 19921) as follows: 



C'rad + Cvd 



)G 



T 



lOiScm-s; V5,000K 



(18) 



where nni and Ue are the number density of neutral hydro- 
gen atoms and that of free electrons. Also, Ciad, C^dw and 
Cstark are natural broadening, van der Waals broadening, and 
linear Stark broadening constants respectively. We adopt this 
para metrisation along with t he va lues of broadening constants 
from ILuttermoser & Johnsonl ( Il992h . For example, Crad = 8.2 x 
10"^ A, Cvdw = 5.5 X 10"^ A and Cstark = 1-47 x 10"^ A for 
Hq. 

To find an observed line profile, we compute the flux as a func- 
tion of frequency. The integration of the flux is performed in the 
cylindrical coordinate system (p, q, t) which is obtained by rotat- 
ing the original stellar coordinate system (p, (j), z) so that the t-axis 
coincides with the line of sight of an observer. The observed flux 
{Fv) is then given by: 



p sin qlvdqc 



(19) 



where d, Pmax, and F are the distance to an observer, the maxi- 
mum extent to the model space in the projected (rotated) plane, and 
the specific intensity {F) in the direction on observer at the outer 
boundaiy. For a given ray along t or alternatively using the optical 
depth {tv defined later in equation l2ll as an integration variable, 
the specific intensity is given by: 



F = he. + j S„ (tI) e ^" drl 



(20) 



where Jo and Sv are the intensity at the boundary and the source 
function (equation |16t of the stellar atmosphere/wind at a fre- 
quency V. The upper integration limit r^a is the total optical depth 
between an observer and the initial integration point. For a ray 
which intersects with the stellar core, Jq is computed from the 
stellar atmosphere model of lKuruc3 ( 1 19791) (cf. Section B. l.lt . and 
7o = otherwise. If the ray intersects with the hot ring on the stel- 
lar surface created by the accretion stream, we set /o = By (Tring) 
where is the Planck function and Tring is the temperature of 
the hot ring. The initial position of each ray is assigned to be at the 
centre of the surface element {dA = p sin q dq dp). We numeri- 
cally integrate equation ( |19t . A linearly spaced radial grid is used 
for the area where the ray intersects with magnetosphere, and a log- 
arithmically spaced grid is used for the wind and the accretion disc 
regions. The numbers of grid points used for a typical integration 
are Up = 180, Uq = 100, and n„ — 101. The optical depth in 
equation[20]is defined as: 



X. (t') dt' 



(21) 



where Xv is the opacity of media, and toe is the total distance be- 
tween the initial integration point and the observer (or to the outer 
boundaiy closer to the observer). Equation (I2U is also numerically 
integrated by using the intersection points of a line along t and the 
AMR grid cells, in which opacities values (xi/) are stored. For a 
high optical depth, additional points are inserted between the orig- 



inal points along the line, and emissivity and opacity are interpo- 
lated to those points to ensure dr^ < 0.05 for the all line segments. 



3 CODE TESTS 

3.1 ID test : comparison with cmfgen model of 
i Hillier & Miller (1998) 

We first test the validity of our statistical equilibrium calcula- 
tion and the observed profile calculations in a simple geometry. 
For this puipose, we consider a case of spherical stellar wind 
from an Of supergiant star since it often exhibits hydrogen and 
helium in emission. A detail atmosphere model of AV 83 (an 
07 laf-l- supergiant ) in the Sma l l Mag ellanic Cloud (SMC) has 
been presented by Hillier et al.l (l2003h . Their models are com- 
puted by the non-LTE line-blanketed atmosphere code CMFGEN 
failher & Milleilll998h in which the level populations are solved 
in the co-moving frame and in a spherical geometry. This code 
has been well-tested and us ed in numerous investigations of stellar 
winds of massive stars (e.g. Hillier & Miller 1999; Crowther et aU 
1999l:lDe Marco et al.! 2000: Dess art et al..2000; Herald et al.,20o'il : 
Mar tins. Schaerer. & H illier 20051: iHillier et al.l2006h . AV 83 is one 



of t wo O stars in the SMC which shows a n ormal Of spectrum 
(cf. IWalbometal]|l995l : IWalbom et al1l2000h . It has a relatively 
dense stellar wind, and consequently exhibits m any w i nd fe atures 
in their spectrum. The model presented bv Hillier et al. 1 (120031) sug- 
gests that AV 83 has the effective temperature of ~32,800 K and 
the mass-loss rate of 7.3 x 10~^ Mq yr~^. 

In the following, we compare the results of the statistical 
equilibrium calculations and the line profiles computed by TORUS 
and CMFGEN. Both models assume the hydrogen and helium 
abundances X = 0.834 and Y — 0.167 (by numbers), re- 
spectively. To match the same physical conditions in the stellar 
wind structure, TORUS adopts the temperature, density and ve- 
locity of the wind used in CMFGEN, which are shown in Fig. |2] 
We also adopt the continuum flux from CMFGEN in comput- 
ing the radiation fields in equations [TT] and [12] Further, we con- 
sider the case without the wind clumping, i.e. a smooth stel- 
lar wind case. The CMFGEN model data of AV 83 used here 
for comparison are available on the web site with the address: 
,http:/ /kookaburra.phyast.pitt.edu/hillier/web/CMFGEN.htm 

Since our code TORUS has newly adopted helium ions in the 
statistical equilibrium calculation, here we concentrate on the com- 
parisons of the helium calculations. The hydrogen level popula- 
tion calculations had been already test ed elsewhere ( lHarriesiri995l : 
lHarriesll2000t ISvmington et al.ll2005l) . Fig. |2] shows the number 
fractions of helium ions in three stages (He I, He II and He III) as a 
function of radial distance from the centre of the star (r'). The dom- 
inating ions of helium are He II and He III, and the number fraction 
of He I typically ~ 10* to ~ 10^ times smaller than those of other 
ions. The figure shows the agreement of the ion fractions computed 
by both codes is good at all radii, except at r' ~ 2 (measured in 
stellar radius) where He I and He II fractions of TORUS are smaller 
than those of CMFGEN by a factor of a few. 

As one can see in Fig. |2l in this particular model of AV 83, 
He II is the most abundant ion up to a several radii. Correspondingly 
the optical thickness of the continuum photon beyond the He II Ly- 
man continuum is expected to be very high. This contradicts with 
our assumption made in Section |Z2l which stated that the gas is 
optically thin to the photospheric continuum radiation. Therefore, 
when the He II Lyman continuum is optically thick, we introduce 
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Figure 2. Left panel: The density (top), temperature (middle) and radial velocity (bottom) of the hot stellar wind model of AV 83 computed by CMFGEN. The 
values are shown as a function of radius r' , which is in the unit of stellar radius (19.6 Rq). The same density, temperature and velocity distributions are used 
in TORUS for the calculation of the statistical equilibrium. Right panel: Comparison of the number fractions of helium ions as a function of radius, computed 
by CMFGEN (solid) and TORUS (dashed). The fractions here are defined as Afjon/Aftotal where A^ion is the number density of helium ions (He I, He II or 
He III), and A'total is the total number density of ions, including both hydrogen and helium. The fraction of He I is much smaller compared to those of He II 
and He III. Overall agreements of the two codes are very good. 



the following special treatment in the statistical equilibrium. We 
assume that the radiation field beyond the He II Lyman edge will 
be reduced to the local source function, and there is radiative detail 
balance in the He II Lyman continuum. In practice, we simply set 
the photorecombination and photoionization teiTns in equations |4] 
and [5] for the ground state of He I to negligibly small values. The 
atmosphere is still assumed to be optically thin to continuum in all 
other wavelengths. The same special treatment of the He II Lyman 
continuum was used by KC78 who modelled the winds from Of 
stars. Without this treatment, our ionization fraction of helium ions 
will be very different from those of the CMFGEN model at small 
radii (e.g. r' < 2) . 

Fig-EJcompares the line emissivity and opacity of He I A10830 
computed by TORUS and CMFGEN as a function of radius. As men- 
tioned earlier, this particular line has been identi fied as one of 
the important wind diag nostic lines in CTTSs (e.g. lEdwards et afl 
llOOd : [Kwan et alJllOOl) . The figure also shows the electron scat- 
tering opacity, and the continuum opacity and emissivity at the line 
centre of He I A10830. Overall agreement between the two codes is 
excellent for the continuum emissivity and opacities at all radii. The 
agreements are also good in the line opacity and emissivity except 
for the line emissivity (?7l) at the inner most radii (r' < 1.1) where 
the outflow speed is subsonic; hence, the Sobolev approximation is 
not suitable there. The figure also shows the corresponding optical 
depth (tl) evaluated at the line centre of He I A 10830. 

Using the opacities and emissivities above, we now compute 
the line profiles as seen by an observer. Fig.|4]shows the normalized 
profile of He I A10830 and Ha computed by TORUS and CMFGEN. 
For a simplicity, the broadening effects such as the Stark broad- 
ening and rotational broadening are not included in these calcula- 
tions. He I A10830 shows a narrow emission near the line centre 
and a shallow but very wide wind absorption in the blue side. On 
the other hand. Ha profiles show a wider emission and a relatively 
narrow wind absorption in the blue. The extent of the red wings in 
the He I A10830 and Ha lines is similar, and they both reach the 




r' r' 

Figure 3. Comparison of opacities, emissivities and optical depth as a 
function of radius, computed by CMFGEN (solid) and TORUS (dashed). The 
line opacity (xl) ™d emissivity (rjx,) are those of He I A 10830. The con- 
tinuum emissivity (rjc) and the thermal opacity (xth) are evaluated at the 
line frequency of He I A10830. The electron scatteiing opacity is shown as 
Xes. The opacities and emissivities shown are in arbitrary units. The optical 
depth (Ti) is evaluated at the line centre of He I A 10830, and it is the sum 
of the hne and continuum optical depths. The radius (r') is in the unit of 
stellar radius which is 19.6 Rq . The agreements between the two codes are 
excellent for the continuum emissivity and opacities at all radii. The agree- 
ments are also good in the line opacity and emissivity except for those at 
the inner most radii (r' < 1.1) where the optical depth is very high and the 
outflow speed is subsonic. 
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Figure 4. Comparison of line profiles (He 1 A10830: right panel and Ha: 
left panel) computed by CMFGEN (soUd) and TORUS (dashed). The line 
strengths of the profiles computed by TORUS are slightly stronger than those 
of CMFGEN; however, overall agreement of the models is excellent. In both 
lines, the broadening mechanisms such as the Stark broadening and the ro- 
tational broadening are not included. 



velocity ~ 800 km s^^ which is very similar to the terminal veloc- 
ity of the stellar wind (900 km s^^) used in the models. The figure 
shows the emission components of both Ha and He I A 10830 pro- 
files computed by TORtJS are slightly stronger than those by CMF- 
GEN; however, overall agreement of the line profiles with those of 
CMFGEN is very good. 

3.2 2D tes t: comparison with models from lMuzerolle et"ai] 
<200lh 

Next, we examine a slightly more complex geometry, namely the 
axisymmetric magnetospheric accre tion flows on to the CTTS as 
described by Hartm ann et alj ( Il994l) . In this model, the gas accre- 
tion on to the stellar surface from the innermost part of the accretion 
disc occurs through a bipolar stellar magnetic field. The magnetic 
field and the gas stream lines are assumed to have the following 
simple form 

r = rm sin^ 6 (22) 

(cf. iGhosh. Pethick. & LambI [19771) where r and 6 are the radial 
and the polar components of a position vector in the accretion 
stream. The symbol rm is the radial distance to the field line on 
the equatorial plane (6 = 7r/2), an d its value is restricte d to be 
between rin and rout (cf. Fig. 1 in lHartmann et al.lll994t) which 
control the size and width of the magnetosphere. Using the field 
geometry above and conservation of energy, the velocity and the 
density of the a c cretin g gas along the stream line are found as in 
iHartmarm et al.l ( Il994l) . We al so adopt the temperat ure structure 
along the stream lines used bv IHartmarm et al] ( 1 19941) . In fact, the 
flow geometry and the density structure described here are used to 
construct the AMR grid cells shown in Fig. [T] (Section 12. U . The 
same magnetosph eric accretion structure is used in the radiative 
tran sfer models bylMuzerolle et all (l200ll) . ISvrnington et all j2005l) 
and Kurosawa et al. ('200^. 

Muzerolle et al. (2001) presented hydrogen line profiles of 
CTTS with many different combinations of magnetosphere ge- 
ometries, flow temperatures and mass-acc retion rates. As our tes t, 
we consider one representative model of iMuzerolle et al] (1200 ih . 
and compute line profiles to check if we can reproduce their re- 
sults. The following model configuration is used in this test. The 
mass (M,) and radius (i?*) of the central star are assumed to 
be 0.5 Mq and 2 Rq, respectively. The effective temperature of 
the star is set to 4,000 K, and the maximum t emperature in the 
magnetosphere (Tmax in iHartmann et alj|l994h is 7,000 K. The 



Figure 5. Comparison of Ha (left) and Pa /3 (right) profile s for a n ax- 
isymmetiic accretion flow as described by iHartmann et al ] HHi). The 
profiles compu ted by TORUS (solid line) are compared with those from 
iMuzeroUe etalj<200ll) (circles). The mass and radius of the central star are 
= 0.5 Mq and = 2Rq, respectively. The effective temperature 
of the star is 4,000 K. The accretion stream has and the maximum temper- 
ature Tniax=7,000K, and the mass-accretion rate Ma = 10~^ Mq yr~^. 
The accretion stream starts in the disc plane between the radial distances 
rin = 2.2 R, and rin = 3.0 iJ» . The star is viewed at the inclination 
angle of 60°. 



geometry of the magnetosphere is determined by setting rm to 
Tin = 2.2 i?* and rout = 3.0 R, in equation | 22 | Thi s corresponds 



to the "small/wide" model in IMuzerolle et alj l l200lh . Finally, the 
mass-accretion rate is set to 10~^ Mq yr~^. 

Fig. m shows Ha and Pa/3 profiles computed with this 
model setup. The fig ure also shows the profile computed by 
IMuzerolle etal] ( 1200 1[) for comparison. In computing the Hq pro- 
file, we have taken into account the line broadening effect as de- 
scribed in section 12.31 Although the line wings of our Hq pro- 
files are slightly wider than those of Muzerolle et al. (2001J, over- 
all agreement between our profiles and t heir model profiles is ve ry 
good. A similar test can be also found in lSvmington et all ( l2005h . 

We are n ot able to compare our model He I line profiles with 
the models of Muzer olle et ai] ( I2OOIII since their model does not 
include helium atoms. In fact. He I lines would not be formed with 
this model configuration because of the low temperature in the ac- 
cretion stream and/or the lack of high energy (photoionization) ra- 
diation source which could excite or ionize helium atoms. The im- 
portance of the photoionization by high energy photons will be dis- 
cussed further in Sections |4T3]|4.2.2| and |4.3.2| 



4 APPLICATIONS TO CLASSICAL T TAURI STARS 

We now apply the radiative transfer model, described and tested 
earlier, to more comprehensive flow configurations of CTTSs 
which include both a magnetospheric accretion and a wind. In the 
following, we briefly describe our model configurations including 
the sources of the continuum radiation, the disc wind model and 
the stellar wind model. The example profiles are computed with 
the combinations of (1) magnetospheric accretion and disc wind, 
and (2) magnetospheric accretion and the stellar wind. Here, we 
consider only the axisymmetric cases. In both cases (1 ) and ( 2), 
the magnetospheric accretion model of IHartmann et al.l d 19941) is 
adopted, which is briefly described in Section |T2l and its accretion 
geometry is shown in Fig.[T] 
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4.1 Continuum sources 

4.1.1 Photosphere 

We adopt stellar parameters of a typical classical T Tauri star for 
the central continuum source, i.e. its stellar radius 7?* — 2.0 Rq 
and its mass M, =0.5 Mq. Consequently, we adopt the effective 
temperature of photosphere Tph = 4, 000 K and the surface gravity 
log Qt = 3.5 (cgs), and use the model atmosphere of Kurucai 19791) 
as the photospheric contribution to the continuum flux. 



4.1.2 Hot spots/rings 

Additional continuum sources to be included are the hot spots/rings 
formed by the infalling gas along the magnetic field on to the stellar 
surface. As the gas approaches the surface, it decelerates in a strong 
shock, and is heated to ~ 10^ K. The X-ray radiation produced 
in the shock will be ab sorbed by the gas locally , and re-emitted 
as op tical and UV light jCalvet & Gullbring|ll998l : lGullbring et al.l 
l200d ) - forming the high temperature regions on the stellar sur- 
face with which the magnetic field intersects. In this study, we 
adopt a single temperature model for hot spots/rings by assum- 
ing the free-falling kinetic energy is thermalised in the radiating 
layer, and is re - emitte d as blackbody radiation, as described by 
iHartmann etal.l ( Il994l). The same model is used in the radiative 
tran sfer models bylMuzerolle et al.l (l200ll) . ISvmington et al] ( l2005l) 
and lKurosawa et al] ( l2006h . The hot spot/ring temperature {T]^s) in 
this model is summarised as: 



4:71 Ria (cos 6in — cos ^ou 



(23) 



where Lhs and a are the hot spot/ring luminosity and the Stefan- 
Boltzmann constant. The angles 6lin and Sout represent the posi- 
tions (the polar angles) where the innermost and outermost dipolar 
magnetic field line intersects with the stellar surface (cf. Figs[T] 
|6l |9). Using equation l l22t along with the inner and outer radii of 
the magnetosphere (-Rmi and -Rmo, respectively), we find 6lin = 
sin-i [iR,/R^i)]^^^ and 6lout = sin"! [(7?,/i?,^o)]''''. The hot 
spot/ring luminosity in this model can be written as 



rI 



1 - 



2R^ 



{Ryni + R n 



(24) 



where Ma is the ma ss-accretion rate. Readers are referred to 
IHartmann et al.l ( Il994h for details. With the stellar parameters used 
in Section l4m along with Ma = 10"^ Mq yr"\ i?nii = 2.2 R, 
and Rmo = 3.0 R,, the corresponding hot spot/ring temperature 
and the accretion luminosity are Ths = 6, 400 K and Lhs = 
1.9 X lO^^ergs"^ = 0.49 Lq. 

Alternatively, we can also construct the multi-temperature 
hot spot mode l obta ined in the 3D MHD simulations by 
iRomanova et al. I j2004l) in which the local hot spot luminosity is 
computed directly from the inflowing mass flux on the surface of 
the star. In this model, the temperature of the hot spots is deter- 
mined by conversion of kinetic plus internal energy of infalling 
gas to a blackbody radiation. They obtained the position-dependent 
mass flux crossing the inner boundary; hence, achieving a position- 
dependent temperature of the hot spot, which can be written as 



1 2 
-V 

a \ 2 



1/4 



(25) 



where p, Vr and v are the density, the radial component of veloc- 
ity and the speed of gas/plasma, respectively. The specific enthalpy 



of the gas is w = 7 (p/p) (7 — 1) where 7 is the adiabatic in- 
dex and p is the gas pressure. This model is useful when the ac- 
cretion flow from a MHD simulati on (e.g. Romanova et alj |2003l : 
iRomanova^et al. 2009) is used, as in lKurosa wa et al. (2008). 

We compare this temperature Ths with the effective temper- 
ature of photosphere Tph to determine the shape and the size of 
hot spots. When Ths > Tph, the location on the stellar surface 
is flagged as hot. For the hot surface, the total continuum flux is 
the sum of the blackbody radiation with Ths and the flux from the 
model photosphere mentioned above. The contribution from the in- 
flow gas is ignored when Ths < Tph. 



4.1.3 X- ray em ission 

CTTSs show relatively strong X-ray emission with their lumi- 
nosities Lx (0.3 — lOkeV) ranging from ^ lO^^ergs"^ to ~ 
10^^ ergs"^ (e.g iTelleschi et all I200I lOudel & Telleschill2007l: 



lOiidel et al.l I2OI0I) . The soft X-ray emission component (with 
the electron temperature To ^a few 10^ K) are possibly as- 
sociated with the shoc k regions either ne ar the base of the ac- 
cretion columns (e . g. ICalvet & GuUbring Il998t iLamzinl Il998t 
iKastner et al] I2OO2': 'Orlan do et aL.20ia) or the base of the stel- 
lar wind/jet (e.g. Gudel et a l.l2008l : ISchneider & Schmitj2008l) . On 
the other hand, the higher temperature (Tc ~ 10 x 10® K) emis- 
sion may be associated with th e flare s and active coronal region^ 
(e.g. Feigel; 



slson & Montm eri3 1 19991: [Ge 
isch et al. 2005; Jard ine et al.l2C 



iGa gne. Skinner, & Danie. 
l2006i;,Stassun et al..200i 



2004'; Preibisch 
Argiroffi, Maggio, & Peres IQQl). 

The coronal heating and the production of X-ray are closely 
related to the for mation and the acceleration of the stellar wind 
from CTTSs (e.g. ICranmeill2009l) . Exploring possible formation 
mechanisms of the stellar wind is important because its mass-loss 
rate is directly related to the amount of the angular momentum that 
can be removed fro m the accreting CTTSs (e.g. iMatt & Pudiitj 
I2005l ; l2007al ; l2008ah . The X-ray irradiation on the accretion disc 
is also important for its thermodynamics, structure and chem- 
istry, and it is the main source of ionization in the outer layers 
of the disc (e.g. Glassgold, Najita, & Igea 2004; Alexander 200^; 
iBergirJ 120091 ; see also lShang et al.ll2002l for the X-ray heating of 
the X-wind). It can also assist a disc to lose its mass at large 
radii via photoevapolation (e.g. lErcolano, Clarke, & Drakell2009l ; 
lOwen. Ercolano. & Clarkell2010l) . The importance of the X-ray ra- 
diation in determining the physical conditions of the circumstel- 
lar material of CTTSs is clear; however, so far, few spectroscopic 
models of optical and near-infrared He and H lines that include 
the effect of photoionization due to X-ray have been explored. 
For this reason, we will include the effect of the X-ray radiation 
in our calculations. No te that a recent local excitation model of 
iKwan & FischeJ JioT l|) demonstrated that the photoionization of 
He due to 'UV radiation' might be also important for ionizing He. 
However, we concentrate on X-ray radiation here, and leave the in- 
vestigation of the relative importance of X-ray and UV radiation as 
our future study. 

Another important motivation for adding the X-ray source in 
our model comes from our initial investigation of the simultaneous 
modelling of H and He emission lines in the CTTSs. Although not 
shown here, we have found that the temperature of the flows (both 
accretion and wind) must be relatively low e.g. T <10,000 B^^for a 
system with the mass-accretion rate, and it could be slightly higher 

See also Fig. 16 in iMuzeroUe etal] i200ll) for the possible temperature 
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for a lower mass-accretion rate. The line strengths of hydrogen 
lines (e.g. Ha) will be unrealistically strong if the temperature is 
much above this temperature. However, at this relatively low tem- 
perature, the coUisional rates of He I are low, and the normal con- 
tinuum sources mentioned earlier (Sections 14. 1 . 1 I and l4. 1.21 do not 
provide enough radiation to ionize/excite He I significantly. This re- 
sults in no significant emission or absorpt ion in He I A10830, and is 
inconsistent with the observations (e.g. Edwards et al. 112006'). One 
way to overcome this problem is to introduce an extra source of He 
ionizing radiation, e.g. an X-ray emitting source. The X-ray radia- 
tion will photoionize He I and keeps its excited states populated. 

Here, we take a simplest approach in implementing the X-ray 
radiation in our model. Firstly, we assume that the X-ray radiation 
arises uniformly from the stellar surface as though it were formed in 
the chromosphere. This assumption allows us to include the X-ray 
emission by simply adding it to the normal stellar continuum flux, 
discussed in Section 14.1.1 1 The dependency of the line formation 
on the different X-ray emission locations is beyond the scope of 
this paper, but will be explored in the future. Secondly, we assume 
that the energy distribution is flat in the energy range between 0. 1 
to 10 keV, and the X-ray flux is normalized according to the total 
X-ray luminosity Lx (0.1 — 10 keV), which we set as an input pa- 
rameter. Alternatively, we could adopt a more realistic energy dis- 
tribution simil^^tojhose used in the X-ray photoevapolation disc 
model of lErcolano et al.|j2008h , who used the synthetic X-ray spec- 
tra of a typical CTTS by assuming the coronal X-ray emission is 
produced in a collision dominated optically thin plasma. The atten- 
uation of X-ray flux is not taken into account in our model; how- 
ever, this can be improved in a future study. A further discussion of 
the X-ray flux attenuation is given in Section [5!2] 

Using the X-ray emission as assumed above, we find a typi- 
cal value of the He I photoionization rate (7Hei) from the ground 
state, at a typical distance from the star (r = 4.0 7?«), is 7hci ~ 
3 X 10~^ when the X-ray luminosity is assumed to be Lx ~ 
2 X lO^^ergs"^. Interestingly, this is v ery similar to the value s 
used in the local excitation calculations o f lKwan&Fis"che3j201lh . 
i.e. 7Hei = 10"" and 10"^ s-\ 

4.2 Disc wind model 

4.2.1 Configuration 

The disc wind mode l used in this analysis is essentially same 
as in [Kurosawa et alj (l2006h. who adopted a simp le kinematical 
wind model of iKnigge. Woods. & Drewl Il995l) and follows 
the basic ideas of the magn eto-centrifugal wind paradigm 



broadly represent the MHD wind 


models (e.g. 


Shu et all 


19941; lUstvueova et al.1 1 19951; lOuved & Pudritz 


19971; 


Romano va et al. 


I997I; Ustvusova et al 


I999I; Konisl & Pudritzl 


2OOOI; KrasnoDolskv et alj |2003|; 


Romanova et alj 20091; 


Murohv. Ferreira. & Zannil l2010h. In the following, we brieflv 



describe our disc wind mode l which is an a daptat ion of the "split- 
monopole" wind model by iKnigge et alj jl995l) . In this model, 
the outflow arises from the surface of the rotating accretion disc, 
and has a biconical geometry. The specific angular momentum is 
assumed to be conserved along a stream line, and the poloidal ve- 
locity component is assumed to be simply a radial from "sources" 
vertically displaced from the central star Readers are referred to 



ranges of the axisymmetric magnetcspheric accretion flow, based on their 
hydrogen emission line profile models. 



I Knigge etal] Il995l) and iLong & Kniggd f20o3) for details. See 
Alen car et al.l ( l2005h and lLima et al.l ( 2010l) for an alternative disc 
wind model. 

There are four basic parameters: (1) the mass-loss rate, (2) the 
degree of the wind coUimation, (3) the velocity gradient, and (4) the 
wind temperature. The basic configuration of the disc-wind model 
is shown in Fig.|6] The disc wind originates from the disc surface, 
but the "source" point (S), from which the stream lines diverge, 
are placed at distance d above and below the centre of the star. 
The angle of the mass-loss launching from the disc is controlled by 
changing the value of d. The mass-loss launching occurs between 
-Rwi and _Rwo where the former is set to be equal to the outer radius 
of the closed magnetosphere (-Rmo) and the latter is set to 0.5 au. 

The local mass-loss rate per unit area (m) is assumed to be 
proportional to the mid-plane temperature of the disc, and is a func- 
tion of the cylindrical radius w = (x^ + y^) i.e. 

m{w) (X [Tiw)]-^" . (26) 

The mid-plane temperature of the disc is assumed to be expressed 
as a power-law in w; thus, T oc lu'. Using this in the relation above, 
one finds 



m (w) oc w'' 



(27) 



where p — 4a x q. The index of the mid-plane temperature 
power law is adopte d from the dust radiative transfer model of 
IWhitnevetaLl ( l2003h who found the innermost part of the accretion 
disc has q = —1.15. To be consisten t with the collimated disc-wind 
model of Krasnopolskv et al] ( l2003h who used p = — 7/2, the value 
of a is set to 0.76. The constant of proportionality in equation |27] 
is found by integrating rh from R„i to J?wo, and normalising the 
value to the total mass-loss rate A/dw 

The azimuthal/rotational component of the wind velocity 
W0 {w, z) is computed from the Keplerian rotational velocity at the 
emerging point of the stream line i.e. {wi, 0) — (GM^/wi)^''^ 
where Wi is the distance from the rotational axis (2) to the emerging 
point on the disc, and by assuming the conservation of the specific 
angular momentum along a stream line: 



Va,{w,z) ^VA,{Wi,0) (—) 
\ w / 



(28) 



Bas ed on the canonical P velo city law of hot stellar winds 
(c.f. lCastor. Abbott. & Kleinlll975l) . the poloidal component of the 
wind velocity (Vp) parametrised as: 



Vp {Wi,l) = Cs {wi) + [f Vcsc — Cs (wi)] ( 1 - 



Rs 



l + Rs 



(29) 



where Cg, /, and I are the sound speed at the wind launching point 
on the disc, the constant scale factor of the asymptotic terminal 
velocity to the local escape velocity (t^csc from the wind emerging 
point on the disc), and the distance from the disc surface along 
stream lines respectively. Rs is the wind scale lengt h, and its value 
is set to 10 -Rwi by following Long & Knig"g3 12OO2I) . We set / = 2 
in the models presented in this paper. 

Assuming mass-flux conservation and using the velocity field 
defined above, the disc wind density as a function of Wi and / can 
be written as 



p{Wi,l) 



m {wi) 



Vp{wi,l) jcos(5| \D{wi,l) cos (5 



d 



(30) 



where D and 5 are the distance from the source point (S) to a point 
along the stream line and the angle between the stream line and the 
disc normal, respectively. The dependency of the density and the 
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Figure 6. Basic model configuration of the disc-wind-magnetosphere liy- 
brid model. The system consists of four components: (1) the continuum 
source located at the origin of the cartesian coordinates (x, y, z) - the y- 
axis is into the paper, (2) the magnetospheric accretion flow, (3) the disc 
wind, and (4) the optically thick but geometrically thin accretion disc. The 
disc wind originates from the disc surface between Wi = R^i and Rwo 
where Wi is the distance from the z axis on the equatorial plane. The wind 
source points (S), from which the stream lines diverge, are placed at dis- 
tance d above and below the star. The degree of wind collimation is con- 
trolled by changing the value of d. The figure is not to scale. 

poloidal velocity on the wind ac celeration parameter [3 is shown in 
Fig. 2 of Kurosawa et alj /2006). Ahhough the outer radius (i?wo) 
of the disc wind could extend to a large radius, most of the mass- 
loss occurs near the inner radius (i?wi) of the disc wind because of 
the rather steep radial dependency of the local mass-loss rate (p = 
— 7/2 in equation[27t adopted here. This concentration of the mass- 
loss near the disc-magnetosphere in teraction region re sembles that 
of the conical wind model of Romanova et all | |2009|) . Further, by 
restricting the outer disc radius (f?wo) to a smaller radius, our disc 
wind model can mimic the conical wind model even better. Lastly, 
the temperature of the disc wind (Tdw) is assumed isothermal as in 
iKurosawa et all bood) . 

4.2.2 Example profiles 

Fig. [7] shows the example helium and hydrogen profiles computed 
for the standard axisymmetric magnetospheric accretion funnels 
(7?mi = 2.2 i?,, Rrno = 3.0 R*, Ma = IQ-^MQyr"^ and 
Tmax = 7, 500 K) and the disc wind with the following param- 
eters: Tdw = 9,000K, Mdw = lO-*M0yr-\ /3 = 0.5, 
i?wi = 3.0 i?., Tiwo = 54,R*,d = 14 7?. (cf. Fig.^. The pa- 
rameters for the central star are as described in Section 14.1 .11 The 
effective temperature of the hotspot, as described in Section |4.1.2| 
is approximately 6,400 K. The X-ray luminosity of the chromo- 
sphere (cf. Section H. 1.3t is assumed to be Lx = 2 x 10^^ erg s~^. 
The profiles shown here are computed for the system with the incli- 
nation angle of z = 45°. The stellar and disk occultations are taken 
into account in our profile calculations. 



We have adopted a relatively low disc wind temperature 
(9, 000 K) here to avoid unrealistically strong emission in the hy- 
drogen lines, produced at the base of the wind. For the same rea- 
son, the temperature of the magnetosphere cannot be much higher 
than Tmax = 7, 500 K (for the system with Ma = 10"^ Mq yr"^ 
or for the corresponding density in the accretion flow). As men- 
tioned earlier, at these relatively low kinetic temperatures of the 
gas, the coUisional rates of He I are not strong enough to populate 
the excited states; hence, an important wind diagnostic line such as 
He I A 10830 will not be formed in emission nor in absorption with- 
out an additional source of excitation/ionization mechanism. For 
this reason, we have added the X-ray continuum (as described in 
Section 14.1.31 ) which can photoionize He I. We find that in general 
the line strength of He I A10830 is sensitive to the amount of the 
X-ray luminosity (Lx) adopted in our model. The model profiles 
computed at three different inclination angles (i = 20°, 50° and 
80°) are summarised in Fig. [7] 

At an intermediate inclination angle (i = 50°), the 
He I A10830 profile in Fig.|7]shows the absorption features in both 
red and blue wings. The relatively narrow blueshifted absorption 
component is caused by the disc wind (outflow), and the absorp- 
tion in the red wing is caused by the magnetospheric accretion fun- 
nel (inflow). The narrow blueshifted absorption component is com- 
monly fo und in the observed H e I A 10830 (e.g. ^ 30 per cent of the 
sample in lEdwards et al .120061) . Similarly, the redshifted absorption 
component is also commonly (47 per cent) seen in the observed 
Hex A 10830 profiles of Edwards el al, (20i)6}. 

The figure shows that He I A5876 profiles are much weaker 
than He I A10830 profiles, with this particular set of model parame- 
ters. The line ratio of He I A5876 and He! A 10830 (using their peak 
fluxes) is about ^ 0.07 for the line profiles computed at the inter- 
mediate inclination angle i = 50° . This is very s imilar to the line 
ratios found in the low temperature gas models of iKwan & Fische3 
( I2OIII) (e.g. see their Figure 9). The relative strength of He I A5876 
is expected to be greater in a model with higher temperature gas 
and higher density than those used here. Our model profiles show 
the redshifted absorption components for all the inclination angles 
even although the absorption is very weak for a low inclination an- 
gle model. On the other hand, the observations of He I A5876 pro- 
files, from 31 CTTSs obtained by Beristain et al. (2001), show that 
the redshifted absoiption is rather rare (3 out 31). The occurrence 
of the redshifted absorption in the model may be decreased when a 
smaller mass-accretion rate was used (e.g. 10~* Mq yr~^). Both 
our models and the observations do not show the blueshifted ab- 
soiption in He I A5876. 

The Ha profile shows broad wings (extending up to 
±400 km s~^) which are mainly caused by the broadening mech- 
anism explained in Section 12.31 It also shows a relatively nar- 
row blueshifted absorption component, caused by the disc wind. 
However, unlike He I A 10830, the redshifted absorption below the 
continuum is not seen in this Ha profile. Based on the morphol- 
ogy of the profile, the profile is classified as 'Type III-B' accord- 
ing the classification scheme of Reipurth et al. (1996). This type 
of line profile is fou nd in about 33 per cent of the observed Ha 
( lReipurthetai]|l996l) . 

Interestingly, the wind signature (the blueshifted absorp- 
tion component) is absent in the model profiles of Pa/3 and 
Br7. In both lines, only the redshifted absorption coinpo- 
nent, caused by the accreting gas in the magnetosphere, is 
present. This type of line profile is often referred to as the in- 
verse P-Cygni (IPC) profiles, and is well reproduced b y the 
magnetospheric accretion model of iHartmann et al.l ( 1 19941) (see 
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Figure 7. An example model calculation with a combination of axisym- 
metric magnetosphere and disc wind. The normalized model profile are 
shown as a function of velocity. The model parameters adopted here are: 
Tmax =7,500 K, Tdw =9,000 K, Lx = 2 x lO^S ergs-\ and Ma = 
lO~'^M0yr~l, Afdw = lO~*M0yr"i (see text for other parame- 
ters). The profiles are computed at inclination angles i = 20° (dashed), 
50° (solid) and 80° (dotted). N.B. the lower limits of the vertical axis for 
He I (10830) and Ha are set to 0, but those for He I (5876), Pa/3, Pa7 and 
Br7 are set to 0.5. 
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who found an almost complete absence of blueshifted ab- 
sorption components and a relatively high frequency of IPC pro- 
files in Pa/3 and Br7. Similarly, the wind signature is absent in 
Pa7 profiles shown i n Fig. |7] This agrees with the observation of 
lEdwards et al.l ( l2006h who presented the atlas of Pa7 profiles from 
38 CTTSs, which also shows little or no signature of wind absorp- 
tion components. The moiphology of the Pa7 model profiles shown 
here overall ag rees with the ra nges of the profile shapes found in the 
observation of lEdwards et al.l(l2006i) . 

As mentioned earlier, the sensitivity of He I A10830 to the in- 
nermost winds of CTTSs and its usefulness for proving the phys- 
ical conditions of the winds have been we ll dem onstrated in the 
past (e.g. Edwards et al.ll2006l ; iKwan et alj|2007l ; lEdwardsll2009l ; 
iKwan & Fischer .201 1 ). To test whether the magnetosphere-disc 
wind hybrid model can account for the types of He I A10830 pro- 
files seen in the observations, we have run a small set of profile 
models for a various combinations of A/dw, T^w and i. The rest 
of the parameters are kept same as in our canonical model shown 
in Fig.|7l We then examined if there is a resemblance between any 
of our model profiles to a set of He I A 10830 line observations in 
lEdwards et Z| ( |2006|) . In the following, we present an example case 



in which the model profile morphology resembles that of an obser- 
vation. 

Fig. [8] shows a simple comparison of a He I A 10830 line 
profile from our magnetosphere-disc wind hybrid model with the 
ob served He I A 10830 profile of UY Aur (CTTS; MO) obtained 
by lEdwards et al.ll200Q . The model profile in the figure is com- 
puted with the following parameters: Tdw = 9,000 K, Mdw = 
2.5 X 10~®MQyr~^ and i = 48°. The mass-accretion rate for 
UY Aur adopted here (log Ma = —7.0) is s imilar to the observa - 
tionally estimated value of log Ma = — 7. 2 tedwards et Zll2006l) , 
but the adopted wind mass-loss rate is slightly smaller than the 
observed value of log7l/„ = —8.2 (Edwards et al. 2006). The 
model reproduces the basic features of the observed He I A10830 
line profile remarkably well. We find good matches of the model 
and the observation in; (1) the location and the narrowness of the 
blueshifted absorption line (caused by the disc wind), (2) the line 
strength, and (3) the depth and the width of the redshifted absorp- 
tion (caused by the magnetospheric accretion). This clearly demon- 
strates the plausibility of our model (magnetosphere -I- disc wind) in 
explaining the type of line profile exhibited by UY A ur. Our model 
agrees with the earlier finding bv iKwan et al.l ( |2007|) and Edward^ 
(2009|) who suggest that the relatively narrow blueshifted absorp- 
tion seen in the observed He I A10830 line of UY Aur is likely 
caused by a disc wind. 

In addition to He I A 10830, we have also com puted Pa7 pro- 
file for UY Aur since the observation presented bv Edwar ds et al.l 
( 2006) also includes Pa7 (which is obtained simultaneously with 
He I A10830). The result is shown also in Fig. [8] Both model and 
the observation shows no sign of blueshifted wind absorption. A 
small amount of redshifted absorption is seen in both model and 
the observation. The overall agreement between the model and ob- 
served Pa7 profiles is very good. An important difference between 
the model and observation is in the strength of the blue wing emis- 
sion. The observation shows a very notable amount of emission in 
the wing (V < —100 km s^^). However, the model do not show 
this emission. This may suggest that the emission from the wind 
in UY Aur is much stronger than that of our model. To match the 
extended blue wing, we may need to adjust the wind temperature 
to a higher value. A further discussion on the blue wing emission 
will be given in Section lSTI 

Again, our main objective here is to test whether the 
magnetosphere-disc wind hybrid model can roughly account for 
the types of He I A10830 profiles seen in the observations, but is not 
to derive a set of physical parameters that is required to reproduce 
the observed line profile as this requires more careful analysis. In 
other words, the model profile shown here is not a strict fit to the ob- 
servation. However, given the similarity of our model profiles with 
those of UY Aur, we should be able to derive the physical parame- 
ters of the accretion flow and wind around the object by performing 
the fine grid parameter search and performing the fitting to not 
only He I A10830 line but also to other hydrogen and helium lines, 
simultaneously. This is beyond the scope of this paper, but shall be 
left for a future investigation. 



4.3 Stellar wind model 

4. 3. 1 Configuration 

As sh own in the observations of He I A10830 by lEdwards et al.l 
( l2006l) , about 40 per cent of the CTTSs in their sample 
exhibit P-Cygni profiles. The observations indicate a possi- 
ble presence of rapidly expanding stellar winds, in some 
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Figure 8. Example model calculations with a cotTibination of axisymmet- 
ric magnetosphere atid disc witid. The nomialized model profiles (tight 
panels) are compared with the observed line profiles (left panels) of 
UY Aur obtained by Edwards etal. (2006). The upper and lower panels 
show HelA10830 and Pa7 profiles, respectively. The observed profiles 
are obtained simultaneously (see Edwards et al. 20()6 for detail). The in- 
clination angle used for computing the model profiles is i = 48°. The 
model parameters adopted here 500 K, Tdw =9,000 K, 

Lx = 2 X lO^^ergs-i, Ma = 10~'^MQyr-\ and Afdw = 2.5 X 
IQ-^ Mq yr" ^. N.B. the model line profiles shown here are not strict fits 
to the observations, but are simple comparisons of the profile morphology. 



of CTTSs. The stellar wind threaded by the open mag- 
netic field in the polar region is a possible mechanis m for 
angular momentum loss (e.g. [Ha rtmann & MacGreaoH I1982L 
Mestellll9 84': 'Hartmann & StauffeT 1989; Tout & PringlP 'l992l; 
Paatz & Camenzind 1996; Matt & Pudritz 2005; Matt & Pudritz 



2008a; Matt & Pudritz 2008b; Matt et al. 2010). Recent work by 
Cranmer (2008) and Cranmer (2009) have shown a possible mech- 
anism of heating and acceleration of the stellar wind, in the context 
of CTTSs, driven by accretion. 

We replace the disc wind, described in Section l4.2.1l with the 
stellar wind, which arises only in the polar directions as shown 
in Fig. 121 Our stellar wind model consists of naiTow cones with 
the half-opening angle 0sw in the polar directions {±z direction as 
shown in the figure). The wind is assumed to be propagating only 
in the radial direction, and its velocity is described by the classical 
beta- velocity law (cf. lCastor & Lamerslll979t) i.e. 



. (r) = no + (vo 



1 _ ^ 
r 



(31) 



where Voc and vq are the terminal velocity and the velocity of the 
wind at the base (r — Rq). Assuming the mass-loss rate by the 
wind is Ms„ and using the mass-flux conservation in the flows in 
cones, the density psw of the wind can be written as: 



Psw (r) = 



Ms, 



(32) 



Anr^Vr (r) (1 — cos &s„) 

where 6'sw > 0, and p^w becomes that of a spherical wind when 
^sw = 90°. The temperature of the stellar wind (Tsw) here is also 
assumed isothermal, as in the case for the disc wind model (Sec- 
tion |4.2.U . 

When the wind is launched from the stellar surface (i.e. i?o = 
i?,), the half-opening angle 6sv, is restricted to be less than ~ 35° 



to avoid an overlapping of the wind with the accretion funnels (cf. 
Fig- Is}. On the other hand, if we assume that the wind is launched 
at a slightly larger, e.g. at the outer radius of the magnetosphere 
Ro = Rmo, the opening angle of the wind is not restricted by the 
geometry of the accretion funnels. In the models presented in this 
paper, we adopt Rq = i?mo(= 3.0 i?,) and Ssw = 50°. 

4.3.2 Example profiles 

Fig.[TO]shows the example helium and hydrogen profiles computed 
for the standard axisymmetric magnetospheric accretion funnels 
(i?mi = 2.2 R,, Rmo = 3.0 7?*, Ma = lO'^MQyr"^ and 
Tmax ~ 7, 500 K) and the stellar wind with the following pa- 
rameters: Ts„ = 8,000K, Mw = 10"*MQyr-\ 0sw = 50°, 
vo = 10kms"\ t)oo = 400kms"\ Ro = 3.0 i?*, and ^ = 0.5 
(cf. Section [4. 3. II and Fig. |9). The parameters for the central star 
are as described in Section 14.1.11 The effective temperature of the 
hotspot, as described in Section 14.1.21 is approximately 6,400 K. 
The X-ray luminosity of the chromosphere (cf. Section |4.1.3t is as- 
sumed to be Lx = 2 X 10^^ ergs~^. The profiles shown here are 
computed for the system with the inclination angle of i = 53°. 
Note that the wind half opening angle 6'sw here is much wider 
(~ 50°) than the one shown in Fig. |9] since the base of the stel- 
lar wind in this particular model is set to be at the outer radius 
of the magnetosphere (_Ro — i?mo) instead of the stellar surface 
(Ro = R*). The stellar and disk occultations are taken into ac- 
count in our profile calculations. 

As in the disc wind cases (Section [4. 2. 2t . we have adopted a 
relatively low stellar wind temperature (8,000 K) here to avoid un- 
realistically strong emission in the hydrogen lines, produced at the 
base of the wind. Again, for the same reason, the temperature of the 
magnetosphere cannot be much higher than Tmax ~ 7, 500 K (for 
the system with Ma = lO"'^ Mq yr~^). As done for the disc wind 
case, we have added the X-ray flux (as described in Section |4.1.3l l 
in the stellar radiation as an additional source of photoionization. 
This helps to populate the excited states of He I and to produce a 
reasonable line strength in He I A10830. The resulting profiles are 
shown in Fig. 1101 

The He I A 10830 profiles (Fig. \T^, computed at the low 
and intermediate inclination angles (i — 20° and 50°), show 
the absorption features in both red and blue wings. However, 
unlike the disc wind case, the blueshifted absorption is much 
wider and it extends almost to the terminal velocity of the wind, 
Voo = 400 km s~^. This is the classical P-Cygni profile which is 
caused by the expanding stellar wind. As mentione d earlier, about 
40 pe r cent of the CTTSs in the observations of lEdwards et al.l 
(200^) show the P-Cygni profiles in He I A10830. The model 
profile of the He I A10830 also shows the absorption in the red 
wing, which is caused by the magnetospheric accretion funnel (in- 
flow). Again, the redshifted absorption component is also com- 
m only (47 per cent) fo und in the observed He I A10830 profiles 
of lEdwards etai] ( 12006 ). 

The model line profiles of He I A5876 (Fig.fTOt are very sim- 
ilar to those computed with the disc wind (Fig. |7}. Again, there is 
no clear sign of a blueshifted wind absorption in the models. The 
similarity of the profiles here with those of the disk wind model 
(Fig- 13 and the fact that the physical conditions (the density, tem- 
perature and velocity) of the magnetospheres used here are exactly 
the same as in the models with the disc wind (Fig.|7]( suggests that, 
at least with this particular set of model parameters, the emission 
in He I A5876 mainly occurs in the magnetosphere. As in the disc 
wind case. He I A5876 profiles are much weaker than He I A10830 
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profiles. The line ratio of He I A5876 and He I A10830 (using their 
peak fluxes) is about ^ 0.07 for the line profiles computed at the in- 
termediate inclination angle i — 50° . As mentioned earlier, the red- 
shifted absorption in He I A5876 i s rathe r rare (abo ut 10 per cent) in 
the ob servation of lBeristain et al.l ( l200lb (see also lAlencar & Basril 
I2OOOI) . 

The Ha profile here also shows a broad wing in red, ex- 
tending up to +400 km s^^. The blue wing is strongly affected 
by the presence of the wind absorption component for the mod- 
els computed at i = 50° and 80°. Similar to He I A10830, the 
blueshifted absorption in the Ha profile is wide and it extends al- 
most to Vac = 400 km s~^. This profile can be also classified as 
a P-Cygni profile. Unlike He I A 10830, the redshifted absorption 
does not go below the continuum level in this Ha profile. Based 
on the morphology of the profile, it i s classified as 'Type IV-B' 
according the classification scheme o f iReipurth et al] ( Il996[) . This 
type of line profile is relativ ely rare (~7 per cent) in the observed 
Ha profiles of IReipurth et al.l il99(h . 

As we found in the disc wind case (Section |4.2| l, the wind 
signature (the blueshifted absorption component) is absent in the 
model profiles of Pa/3 and Br7. Only the redshifted absorption com- 
ponent cased by the accreting gas in the magnetosphere is present 
(the IPC profiles). This is consistent with the near infrared spec- 
trosco pic observations of 50 T Tauri stars by iFolha & Emerson! 
( l200lh . The model profiles of Pa/? and Br7 shown here are very 
similar to those of the disc wind model (in Fig. [7} because these 
lines are not affected by the presence of the wind, and the inclina- 
tions (i) used in both models are same. The wind signature is also 
absent in Pa7 profiles shown in Fi g.llOl As mentioned e arlier, this is 
consistent with the observation of lEdwards et al. I( l2006h . Again, the 
morphology of the Pa7 model profiles shown here overall agrees 
with the ranges of th e profile shapes found in the observation of 
lEdwards et^ ( l2006l) . 

As in the disc wind case, we now test whether the 
magnetosphere-stellar wind hybrid model can roughly account for 
the types of He I A10830 profiles seen in the observations. We have 
run a small set of profile models for a various combinations of A/sw, 
Tsw, i, and the turbulence velocity in the magnetosphere and the 
wind. Note that no turbulence velocity is included in the previous 
models. The rest of the parameters are kept same as in our canoni- 
cal model shown in Fig. [TO] We then examined if there is a resem- 
blance between an y of our model profile s to a set of He I A10830 
line observations in lEdwards et al.l ( l2006[) . 

Fig. [TT] shows a simple comparison of a He I A10830 line 
profile from our magnetosphere-stellar wind hybrid model with 
the observed He I A 1 0830 profile of AS 353 A (CTTS; K5) from 
lEdwards et al.l ( |2006[) . The model profile in the figure is com- 
puted with the following parameters: Tsw = 8,000 K, A/sw ~ 
2.5 X 10-**MQyr-\ Lx = 1 x 10^^ergs-\ i = 10°, and 
the turbulence velocity of 80 km in both the magnetosphere 
and the wind. The rest of the parameters are same as those used in 
Fig. [To] Note that the mass-accretion rate for AS 353 A adopted 
here (logAfa = —7.0) is much sm aller than the observ ationally 
estimated value of log A/a ~ —5.4 jEdwards et al.l l2006). but the 
adopted wind mass- loss rate is very sma ller to the observed value 
of logA/w = -7. 6 jEdwards et al.l2006l) . 




By following' Kwan et alj i2007l ). to match the location of the 
emission peak in red (~ 200 km s~^), a relatively high turbulence 
velocity (80 km in both the magnetospher e and the wind) is 
added to the model. Note that lKwan et alj j2007h applied a high tur- 
bulence velocity in the wind only, but not in the magnetosphere — 
their model does not include a magnetosphere at all. As we demon- 



Figure 9. Same as in Fig.|6] but for the (bipolar) stellar wind and magne- 
tospheric accretion hybrid model. The system consists of four components: 
(1) the continuum source located at the origin of the cartesian coordinates 
{x,y,z) - the y-axis is into the paper, (2) the magnetospheric accretion 
flow, (3) the stellai' wind, and (4) the optically thick but geometrically thin 
accretion disc. The wind is launched from a sphere with radius Ro, but 
is restricted within the cones with the half opening angle 6sw- We adopt 
Ro = Rmo{= 3.0 Rt) and 6sw = 50° for the models presented in this 
paper. The density distribution is symmetric ai'ound the ^-axis. The figure 
is not to scale. 



strate later in Section|5T|(Fig.[T2), the line 'emission' mainly arises 
from the magnetosphere in the models presented in this paper; 
therefore, it is not possible to match the peak of the observed profile 
without invoking a significant (80 km s~^) turbulence velocity also 
in the magnetospheric accretion funnels. This essentially arbitrary 
line-broadening (over-and-above thermal and pressure broadening) 
may have a physical basis in the magnetically-channelled nature of 
the inflows and outflows, and further radiative transfer models will 
be necessary to investigate whether it would be possible to deter- 
mine whether the different turbulent broadening may apply in both 
the wind and the accretion streams. If one could find a solution in 
which the wind emission dominates, the high turbulence velocity 
in the magnetosphere assumed here might not be necessary. 

The figure shows that the magnetosphere-stellar wind hybrid 
model reproduces the basic features of the observed He I A10830 
line profile relatively well. The model profile shows a clear P-Cygni 
profile, and the overall agreement of the model and the observa- 
tional line profile morphology is very well. Again, this demon- 
strates the plausibility of our model (magnetosphere -1- stellar wind) 
in explaining the type of line profile exhibited b y AS 353 A . Our 
model ag rees with the earlier finding by Kwan etal] ( |2007|) and 
lEdwards! ( l20()9b who suggest that the relatively wide blueshifted 
absorption (the P-Cygni profile) seen in the observed He I A10830 
line of AS 353 A is likely caused by a stellar wind. 

We have also computed Pa7 profi le for AS 353 A , and com- 
pared it with the observed profile from lEdwards et aP (2006). The 
data is acquired simultaneously with the He I A10830 line. The pro- 
files are also shown in Fig.[TT] No clear sign of wind absorption is 
seen in both model and observed Pa7 profiles. Note that the line 
width of the model Pa7 profile in Fig.[TT]is much larger than that 
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Figure 10. An example model calculation with a combination of axisym- 
metric magnetosphere and a stellar wind. The normalized model profile are 
shown as a function of velocity. The model parameters adopted here are: 
Tmax = 7,500 K, Tsw = 8,000 K, Ma = IQ-^MQyr-i, A/^w = 
10~* M0 yr~^, and Lx = 2 X 10^® ergs"^. The profiles are computed 
at inclination angles i = 20° (dashed), 50° (solid) and 80° (dotted). N.B. 
the lower limits of the vertical axis for He I (10830) and Ho are set to 0, but 
those for He I (5876), Pa/3, Pa7 and Br7 are set to 0.5. 



of the Pa7 profile in Fig.[TO](for i — 20° case) because of the high 
turbulence velocity (80 km s~^) in the magnetosphere used in the 
model for AS 353 A. The model shows a weak redshifted absorp- 
tion, but the observation does not. The line strength of the model is 
noticeably weaker than the observed profile. Another important dif- 
ference between the model and observation is in the amount of the 
emission in the wings of the profiles. The observation shows sig- 
nificant emissions in the wing extending up to V ~ ±400 km s~^. 
On the other hand, our model do not show strong wing emissions. 
This difference was also found in the comparison of Pa7 profiles 
in UY Aur (Fig. [8}. To produce the extended wing emission and a 
stronger line centre flux in our model, a higher wind temperature 
may be needed. A further discussion on the extended wing emis- 
sions will be given in Section lSTT] 

Once again, our main objective here is to demonstrate whether 
the magnetosphere-stellar wind hybrid model can roughly account 
for the types of He I A10830 profiles seen in the observations, but is 
not to derive a set of physical parameters that is required to repro- 
duce the observed line profile as this requires more careful analy- 
sis. For example, the mass-accretion rate used here is much smaller 
than the observed value and it should not be considered as a derived 
model parameter. However, the resemblance of the model profiles 
to those of AS 353 A suggests a possibility for deriving the phys- 
ical parameters of the accretion flow and wind around the object 
by performing the fine grid parameter search and performing the 
fitting to a set of (preferably) simultaneously observed multiple 
hydrogen and helium lines (including He I A10830). 
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Figure 11. Example model calculations with a combination of axisym- 
metric magnetosphere and a stellar wind. The normalized model profiles 
(right panels) are com pared with the observ ed line profiles (left panels) of 
AS 353 A obtained bv lEdwards etal J j2006) . The upper and lower panels 
show He I A 10830 and Pa7 p rofiles, respectively. T he observed profiles are 
obtained simultaneously (see lEdwards et al]|2006l for detail). The inclina- 
tion angle used for computing the model profiles is i = 10°. The model 
parameters adopted here are: Tmax =7,500 K, Ts„ =8,000 K, Lx = 1 x 
lO^S ergs-\Afa = lO"'' Mq yr"!, and Msw = 2 X IQ-^MQyr-i. 
N.B. the model fine profiles shown here are not strict fits to the observations, 
but are simple comparisons of the profile morphology. 



5 DISCUSSIONS 

5.1 Some unresolved issues 

While the blue and redshifted absorption components in the model 
line profiles (e.g. in Figs. |7] and 110) could provide us the physi- 
cal conditions and locations of line 'absorption', the locations of 
line 'emission' are not so straightforwardly understood from the 
model profiles. To demonstrate which gas flow component (mag- 
netosphere or wind) contributes to the emission part of the profiles, 
we compare the He I A 10830 line profiles computed using (1) the 
magnetosphere only, (2) the disc/stellar wind only, and (3) the com- 
bination of the magnetosphere and disc/stellar wind. The results are 
placed in Fig. [12] The physical parameters used here are the same 
as those in Figs. [7] and [TO] for the models with the disc wind and 
that with the stellar wind, respectively. As one can see from the 
figure, with these particular sets of model parameters, mainly the 
magnetospheres contribute to the emission. The wind emission is 
relatively low, and the winds mainly contribute to blueshifted ab- 
sorptions. Although not shown here, the magnetosphere is also a 
main emission contributor in other lines shown in Figs. iTlandfTOl 
Since we have not yet surveyed a large parameter space for the 
wind and magnetosphere in our models, we cannot generalise the 
statement above. The relative importance of the wind and magne- 
tospheric emission and the determination of the emission location 
require further investigation which we defer to a future study. 

There are at least two important differences between the 
model profiles, computed with the magnetosphere as a main line 
emitter, and observed profiles. Firstly, the line widths of Pa/3 and 
Br7 in our models (e.g. Figs. [7] and 110) are much narr ower than 
those from the observation by Ipolha & Emers^ (l200lb . The full 
width half-maxima (FWHM) of Pa/3 and Br7 are about 200 km 
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in iFolha & EmersonI ( 1200 ih . but those in our model are ~ 100 
or less. Since the maximum infall velocity in our model is essen- 
tially limited to the escape velocity from the star (~ 300 km s^^ 
with R, = 2Rq and A'/* = 0.5 Mq), it is hard to achieve the 
FWHM values found in the observation without an additional form 
of line broadening (e.g. turbulence broadening). iKurosawa et al.l 
1 I2OO6I) also discuss a similar problem. The emission from a stronger 
wind may also help to make the lines broader, but the stronger wind 
would also cause a blueshifted absorption. 

Secondly, the extent of the blue wing emission in the model 
profiles are much smaller than those found in observations. In the 
comparison of the model profiles with the observation of UY Aur 
(Fig. [8j, we found the blue wings of the observed He I A10830 
and Pa7 extend up to ^ 300 km while our models show the 
blue wings only up to ~ 130 km s~^. Similarly for AS 353 A 
(see Fig. lilt , the blue wing of the observed Pa7 profile extends 
up to ~ 400 km s^^ while that of the model extends only up to 
200 km s^^, even with a large turbulence velocity (80 km s~^) in- 
cluded. As mentioned earlier, in the models considered here, the 
line emission mainly arises in the magnetosphere. Hence, with this 
set up, the models have difficulties in producing the blue wing 
emission extending to such high velocity. This discrepancy may 
suggest that the emission from a fast wind (a disc or stellar wind) 
may be necessary to reproduce the extended blue wing emission 
seen in the emission. Whether one can reproduce the strong blue 
wing emission by a wind without causing a strong blueshifted ab- 
sorption and a (too) strong central emission component is still un- 
known. This issue will be investigated further in a future study. 

As we have briefly mentioned in Section l4.3.2l Ha lines com- 
puted with the stellar wind + magnetospheric accretion (see Fig.llOt 
exhibit strong P-Cygni like blueshifted absorptions ; however, this 
type of profile (Type IV-B in the classification by iReipurth et al.l 
fl99^ is rare in the observations. In our model configuration (with 
a stellar wind), it is difficult not to produce a P-Cygni like absorp- 
tion when we require a P-Cygni like absorption in He I A10830. 
In other words, whenever a P-Cygni like absorption is present in 
He I A10830, it also appears in Ha. One possible way to avoid this 
problem is to use a lower wind temperature which reduces the ion- 
ization level of hydrogen in the wind, but not that of helium since 
the helium ionization level is mainly controlled by the X-ray flux. 
Another possibility is to increase the temperature of the magneto- 
sphere which would increase the emission in the wings of Ha and 
would suppress the strong blue absorption features. This point shall 
also be investigated further in a future study. 



5.2 X-ray optical depths 

When computing the statistical equilibrium of the atomic popula- 
tion levels, we have neglected the attenuation of X-ray emission by 
implicitly assuming the radiation is optically thin. Here we exam- 
ine whether this assumption is reasonable, by computing the X-ray 
optical depth (continuum) using the ionization structure found in 
the disc wind plus magnetospheric accretion model shown in Fig.|7] 
(Section r4.2.2t and also the stellar wind plus magnetospheric accre- 
tion model shown in Fig. 1 101 (Section [4.3. 2t . Note that the optical 
depths computed from these models are not self-consistent since 
the inclusion of the X-ray attenuation should, in principle, alter the 
ionization structures itself, especially for He I. Hence, the following 
optical depth values are merely rough estimates. 

For the photon energy e = 0.1 keV and along a line of sight 
to an observer at i = 50°, the optical depth across the magneto- 
spheric accretion stream is rx ~ 10^, and that across the wind 
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Figure 12. Compaiisons of model He I A10830 line profiles computed with 
different flow elements in the disc wind -I- magnetosphere configuration (left 
panel) and in the stellar wind -I- magnetosphere configuration (right panel). 
The profiles are computed with (1) the magnetosphere only (dashed), (2) the 
disc/stellar wind only (dotted), and (3) the combination of the magneto- 
sphere and disc/stellar wind (solid). The physical parameters used here are 
the same as those in Figs. [T] and [lO] for the models with the disc wind and 
that with the stellar wind, respectively. The inclination angle i = 50° is 
used for all the models. 



is Tx ~ 20 and ~ 40 for the disc wind and stellar winds, re- 
spectively. On the other hand, for e — 1.0 keV, rx ~ 1 for the 
magnetosphere, and rx ~ 0.01 and ~ 0.02 for the disc and stel- 
lar winds respectively. The optical depths are computed with the 
models with the mass-accretion Ma ~ lO^^MQyr^^ and the 
wind mass-loss rate Mw = 10~® Mq yr^^, which are about 10 
times higher than those of a typical CTTS. If more typical values 
are used (i.e. Ma = 10"* Mq yr"^ and A'/„ = lO"'' Mq yr"^), 
the optical depths would be about 10 times lower than the values 
shown above. In our model, the main source of the continuum opac- 
ity around the photon energy 0. 1 keV is the bound-free opacity of 
He I. 

The high optical depths at 0.1 keV indicates that the most 
of the X-ray at this energy could not penetrate the high density 
magnetospheric accretion streams. On the other hand, the optical 
depth across the magnetosphere at 1.0 keV is marginal (rx ~ 1.0), 
and may penetrate the accretion streams, after its intensity being 
reduced by a factor of e^^ ~ 2.7. This suggests that our He I 
line strength may be overestimated for a given amount of X-ray 
luminosity. In reality, the magnetospheric accretion may not be 
completely axisymmetric, and the accretion would occur only in 
two or more funnel streams as shown by 3D MHD simulations 
(e.g. lRomanova et al.l2003l) . This indicates a presence of large gaps 
between the fuimel flows; hence, a significant fraction of the X-ray 
produced near the stellar surface could still reach to the stellar or 
disc wind by escaping though the gaps, and ionize He I in the wind. 
We plan to implement a proper treatment of X-ray attenuation ef- 
fect in a future model. 



6 SUMMARY AND CONCLUSIONS 

The radiative transfer code TORtJS l lHarriesll2000l:l Kurosawa et al.l 
2004; Symington etal 
Kur osawa et aU 12001 



I2OO5I : 
extended 



Kurosawa et al. 



200d ; 



IS extended to include helium atoms 
( 19 levels for He I, 10 levels for He II and the continuum level for 
He III). Only hydrogen atoms (up to 20 levels) are included in the 
previous versions of the code. The addition of he lium atoms in 
our co de is motivated b y the r ecent observations of lEdwards et"al] 
l l2003l) . iDupree et al.] ( l2005h and lEdwards et al.( 1 I2OO6I) who 
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demonstrated a robustness of optically thick He I A10830 as a 
diagnostic tool for probing the innermost wind of CTTSs. 

We have tested our new implementation of the statistical 
equilibrium routine which simultaneously solves H and He level 
populations (Section |3. Il l by comparing the mod el with the non- 
LTE line-blanketed atmosphere code CMFGEN jHillier & Miiierl 
Il998h in which the level populations are solved in a spherical 
geometry (ID). For the case of spherically expanding hot stel- 
lar wind (AV 83), we have found a good agreement between 
TORUS and CMFGEN in computing the number fractions of He 
ions (Fig.|2l(, opacities (Fig. [3}, and observed line profiles (Fig.|4]l. 
The code has been also tested in a 2D (axisymmet ric) geometry 
(Secti on [3.2| >. The emission line profile models of Muzer olle et al.l 
( l200lh . who used the axisym metric magnetospheric accretion 
model of iHartmann et al.l ( Il994) . have been well reproduced with 
our model (Fig.|5}. 

We have also demonstrated the capabilities of our newly 
implemented code in the line formation problems of CTTSs 
(Section |4ll. We have considered two different combinations of 
the inflow and outflow geometries: (1) the magnetospheric ac- 
cretion and the disc wind (Section 14.2) . and (2) the magneto- 
spheric accretion and the bipolar stellar wind (Section |4.3t . In 
both cases, our models are able to produce l i ne profiles similar to 
the observed ones (e.g. Edwards et al. "l994'; ' Reipurth et al. 
lFolha&Emersonl l2001: Edwards et al. 2006), in the morphology 
and the line strengths (Figs. 171 andllOll. In particular, the compari- 
so n of the model He I A 10830 profiles with with the observations 
of lEdwards et"al] ( l2006h (Figs. [8] and [IB suggests that the rela- 
tively deep and wide blueshifted absoiption component found in 
some CTTSs is most likely caused by the bipolar stellar wind while 
the narrower blueshifted absorption is cau sed by the disc wind 
(cf. lKurosawa et alJl2006l : lKwan et alj2007h . 

During our initial investigation of the simultaneous modelling 
of H and He emission lines, we found (in both the disc wind and 
stellar wind cases) the temperature of the wind must be relatively 
low e.g. T ~10,000 K (for Afa = 10"^ Mq yr"^ and it could be 
slightly higher for a lower mass-accretion rate), and it is likely be- 
low 20,000 K. The line strengths of hydrogen lines (e.g. Hq) will 
be unrealistically strong if the temperature is much higher. On the 
other hand, at this relatively low wind temperature (T ~ 10,000 K), 
the collisional rates of He I are too low, and the normal contin- 
uum sources (the photosphere with T^s ^--^4,000 K and the hot 
spot with Ths ~6,400 K or even with Ths ~8,000 K) do not pro- 
vide enough radiation to ionize/excite He I significantly. To popu- 
late the excited states of He I and hence to produce He I A 10830 
with a realistic line strength (c.f. Figs. [8l and II U. one requires a 
high energy radi ation source such as X-ra y radiation. Interestingly, 
a recent work bv lKwan & FischeJ ( 1201 l|) has shown that UV pho- 
toionization is also a probable excitation mechanism for generat- 
ing He I A10830 opacities. From the ob servations of TW Hya and 
RU Lup, Ijohns- Krull&Hercze3j2007h also reached a similar con- 
clusion. Both our model and the work by iKwan & Fischeil JioTlh 
suggest a need for the extra photoionizing radiation sources. The 
relative importance of the UV and X-ray luminosity is not known 
at this moment. 

The following is the list of a several key questions which are 
yet to be answered: Do the stellar wind and disc wind coexist? 
In which physical condition, one dominates the other? What are 
the possible ranges of temperature and density in the winds? Must 
the winds of CTTS be clumpy (inhomogeneous), as suggested by 
iKwan & FischeJ f20Tlh ? What is the relative importance of the UV 
and X-ray? Where exactly is the location of X-ray emission? We 



intent to investigate and provide answers to some or most of the 
questions mentioned above in our subsequent papers in the future. 

Finally, we also plan to apply our radiative transfer model 
to the results of 2D and 3D MHD simulations of the (inner- 
most) conical w i nd an d magnetospheric accretion of CTTSs by 
iRomanova et al. I ( I2OO9I) . A direct comparison of the MHD sim- 
ulations with observations is difficult, but it can be mediated by 
modelling emission line profiles which could constrain some basic 
physical parameters of the flow such as the geometry, temperature 
and mass accretion and outflow rates. In turn, these parameters can 
provide constraints on the amount of angular momentum transfer to 
and out of the stars. By using the different time slices of the MHD 
simulations, we will be also able to address the issues of the time- 
variability in the l ine and continuum by comparing the models with 
observations (e.g. [Kurosawa et alj2008l) . 
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